Phonon Laser Effect and Dicke-Hepp-Lieb Superradiant Phase Transition in Magnetic 
Cantilever Coupled to a Bose Einstein Condensate 
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We propose a possibility of a phonon laser by coupling a Bose-Einstein condensate to a nanome- 
chanical cantilever with a magnetic tip. Due to the magnetic coupling, atomic spin flips induce 
cantilever motion which can be used to produce a phonon laser. The system is described by the 
equivalent of the Jaynes-Cummings Hamiltonian. By controlling the number of atoms and the pop- 
ulation inversion, one can obtain either a continuous wave (cw) or transient lasing. The two-body 
atom-atom interaction is also shown to coherently manipulate the lasing process. We also show that 
in the strong coupling limit, the same system can undergo a Dicke-Hepp-Lieb superradiant phase 
transition. Exotic phase diagrams can be obtained by tuning the two body atom-atom interaction. 
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I. INTRODUCTION 

There have been strong activities recently in quantum 
'optics' with phonons instead of photons, such as experi- 
ments for the study and control of single phonons [l| , or 
theoretical concepts of a phonon laser. 

The phonon analog of the optical laser has been pro- 
posed in numerous physical systems so far. To name 
a few, particular proposals are based on paramagnetic 
ions in a lattice paraelectric crystals isolated 
trapped ions [|| , quantum wells @ , semiconductors , 
nanomechanic systems , nanomagnets , and ultra- 
cold matter [l2|. From the experimental point of view, 
phonon laser action has been demonstrated in cryogenic 
Al 2 3 :Cr 3+ [prf!], A1 2 3 : V 4+ 0, semiconductor su- 
perlattices |17|. harmonically bound magnesium ions flpl 
and very recently in a compound microcavity system [19j . 

Recently the field of cavity optomechanics has become 
an attractive research topic with Bose-Einstein conden- 
sate [2fj| - l3]| and atomic ensembles [32j - [37j . A cavity 
opto-mechanical system, generally consists of an opti- 
cal cavity with one movable end mirror. Such a sys- 
tem is utilized to cool a micromechanical resonator to its 
ground state by the pressure exerted by the cavity light 
field on the movable mirror. The studies on cavity opto- 
mechanics of atoms show that sufficiently strong and co- 
herent coupling would enable studies of atom-oscillator 
entanglement, quantum state transfer, and quantum con- 
trol of mechanical force sensors. Recent experiments have 
shown an impressive level of coherent control over micro- 
and nanomechanical oscillators. Magnetically coupling 
ultracold atoms to mechanical oscillators creates a unique 
setting where coherent quantum control over all degrees 
of freedom can be achieved [H, [3!| . 

In this paper, we propose a phonon laser that operates 
like a two-level optical laser by using a Bose-Einstein 
condensate of 87 Rb atoms magnetically coupled to a 
magnetic cantilever [38[- We also show that for strong 
coupling, the same system can undergo a Dicke-Hepp- 
Lieb type phase transition [4(il - |42j into a phonon super- 



radiance regime. The two-body atom-atom interaction 
thereby plays a crucial role in the dynamics of the sys- 
tem. 



II. PHONON LASING 




Figure 1: Top figure: Schematic view of the coupling of the 
BEC (in red) to the magnetic cantilever (blue). The BEC 
is at a distance yo from the cantilever. The cantilever per- 
forms out-of-plane mechanical oscillations denoted by a(i). 
The oscillatory component of the magnetic field B r (t) cou- 
ples the magnetic cantilever to the atomic spin F. Bottom 
figure: Hyperfme structure of 87 Rb. Transition from the state 
\F = l,rriF = —1) to \F — l,rriF = 0) leads to the emission 
of phonons. 

The phonon lasing device proposed here essentially 
consists of a gas of ultracold 87 Rb atoms at a distance 
lo above a cantilever resonator (38j with a ferromag- 
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netic tip (fig. 1) , which creates a strong magnetic field 



Gr, 



2^0 I Mm I 



Here, p, m is the magnetic moment 



of the ferromagnetic tip and no is the permeability of 
free space. The magnetic cantilever performs out-of- 
plane mechanical oscillations which transduces into an 
oscillatory magnetic field B r (t). The atomic spin F in- 
teracts with B r (t), where fl = [ibQfF is the magnetic 
moment operator. The ground state hyperfine spin lev- 
els \F — l,mp) of 87 Rb are also shown in Fig.l. The 
energy splitting between adjacent mp levels is given by 
the tunable Larmor frequency u)l = HB\9F\Ba/h. The 
tunability of ul allows one to coherently control the de- 
tuning 6 = uj r — lol, where u r is the mechanical fre- 
quency of the cantilever mode. Near resonance [5 w 0), 
atom-phonon coupling leads to spin flips, between the 
ground level |<jr) = \F = 1, mp = 0) and the excited level 
|e) = \F = l,mp = —1). The transition \g) <^ |e) can 
be decoupled from other mp levels [38]. This coupled 
system of BEC and magnetic cantilever (atoms driving 
the cantilever and the cantilever driving the atoms) leads 
to a kind of positive feedback that arises in all laser like 
system. In order to observe phonon lasing, we must also 
introduce a pumping mechanism to compensate for the 
dissipation in the phonon number and the loss of conden- 
sate atoms. This can be achieved by pumping atoms in 
the excited state |e), so as to maintain a steady popula- 
tion inversion. 



Here, ip g and ip e are the ground and excited state wave- 
function of the condensate. Also uii and Vi(r) (i = g,e 
) are the energies and the trapping potentials respec- 
tively for the ground and excited states of the conden- 
sate, g is the atom-phonon coupling constant taken to be 
real, m is the mass of single atom of the condensate and 
a-gg, flee and a ge are the s-wave scattering lengths cor- 
responding to ground — ground, excited — excited and 
ground — excited atomic states respectively. Here we 
have taken a eg 

We now write 14 



(6) 



$ g (r,t) = VNbo(t)e g (r), 



Mr,t) = VNc {t)Ur), 



(7) 



where bo(t) and co(t) are the annihilation operators for 
the ground and excited state atoms, respectively. Here, 
£ g (r) and £ e (r) are the single particle wave functions 
for the ground and excited state respectively satisfy- 
ing the normalization Yli= g e I ^ r l£»( r )| 2 = 1- Ignor- 
ing counter rotating terms, we get the following second- 
quantized Hamiltonian in terms of the normalized oper- 
ators, b — > ^/Nb and c — > VNc , 



A. Hamiltonian 

The coupled dynamics of the magnetic cantilever and 
the BEC can be described by the Hamiltonian, 



H 



Hu> r a?a + huoblbo + hioic^co + %[GdbQc\ + G*a [ cob\^ 

nK eg blb clc , (8) 
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d 3 r £(r) 



-ti 2 V 

2m 



+ Hu g + V g (r) 



Ur), (9) 



Hatom = / d 3 r^(r)[ 



d 3 r^ e (r)[ 



2m 



-h 2 v 2 

2m 



+ Hu) g + V g (r)]ip g (r) 



+ fiw e + V e {r)]i> e (r), (2) 



H phonon — huJ r (l^ CL, 



(3) 



fkj 1 = / d*re e (r) 



2m 



+ huj e + V ge (r) 
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d 6 r \Ur)\\ 



Ur), (10) 



(11) 



Hatom-phonon = + a) [ / ?/>|(r)V>e {r)d 3 r 



+ / 4(r)i> g {r)d 3 r], 



(4) 
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(13) 



(14) 



3 



B. Phonon laser mean field equations 

We now write down the Heisenberg equation of motion 
for the phonon operator a and the atomic operators &o 
and Cq, 



6o f 



a = —iui r a — — — iG* c^b^ 



= l [ UJo bl + aGcl + ^\bo\ 2 bl 



(15) 



Co 



K eg \c \ 2 bl]~ J&J, 



-i{u-iCQ + aGc + -^|co| 2 cJ 



(16) 



Ke 9 \bo\ 2 cl] 



7 

2 C °- 



(17) 



Here, T and 7 are the damping rates for the phonons 
and the atoms, respectively. Now in terms of the po- 
larization p = byCo/N and population inversion An = 
(| Co 1 2 — |£>o| 2 )/-^: the Heisenberg equations of motion can 
be rewritten in the rotating frame of the phonon fre- 
quency u r as, 



a = --a + G*Np, 



7 

iAuip p + GaAn, 



An = J(An eg - An) - 2[G*a f p + Gap*]. 



(18) 



(19) 



(20) 



Here, 
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Also An eq is the equilibrium value of An. 

After factorization, the mean-field steady state solu- 
tions of Eqns. (18)-(20) lead to a critical number of atoms 
N cr required to support a continuous wave (cw) laser, 
N > N cr = T{Au 2 + 7 2 /4)/(|G| 2 An e97 ). For possible 
experimental values mentioned [38j . N cr = 300/An e? . 
For An eg = 0.2, N cr — 1.5 x 10 3 atoms which is a rea- 
sonable number. 



of oscillations, y/\An eq \N G, is larger than the fastest re- 
laxation of the system 7. From this condition, we can 
estimate the minimum number of atoms N t required to 
observe this transient phenomena as Nt — 7 2 / (\An eq \G 2 ) 
« 1.25 x 10 4 . On the other hand if N t > N > N cr , a sin- 
gle pulse can be produced instead of energy oscillations 
between the cantilever and the atoms. 

The two transient outputs mentioned above are shown 
in Fig. 2 by numerical integration of Eqns.(18)-(20). In 
Fig. 2(a) for N = 1 X 10 5 (thin line), coherent energy os- 
cillations are observed while for N = 4 x 10 3 (thick line), 
one single pulse is obtained which takes away a large 
part of the energy at one go. If the number of atoms N 
becomes less than N cr , the single pulse disappears com- 
pletely. The tails of the output pulses decay as Fig. 
2(b) shows the influence of the two-body interaction on 
the transients. Keeping the number of atoms N = 1 X 10 5 
fixed, increasing the value of {K + /A — K eg ) decreases the 
amplitude of the coherent oscillations but increases the 
frequency of oscillations. Atomic two-body interactions 
can be manipulated by Feshbach resonances [45j . 
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C. Transient solutions 

One of the important predictions of the Jaynes- 
Cummings model [44| are coherent population oscilla- 
tions between an oscillator and a (pseudo) spin, i.e. a 
two-level system. Such energy oscillations can also be 
observed in our current system if the effective frequency 



Figure 2: Plot (a): Plot of A(t) = (t)a(t) as a function 
of time for N = 1 x 10 5 (thin line), N = 4 x 10 3 (thick 
line), Plot (b): Plot of A(t) = a'(t)a(t) as a function of time 
for (S± - K eg ) = (thin line), (£± _ K eg ) = 2.2 x 1(T 3 
(thick line). All parameters are dimensionless with respect 
to the atomic damping rate 7. Other parameters used are: 
r/7 = 0.5, G/j = 0.02, u) L - u) r = 0, An eq = 0.2. 
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III. PHONON SUPERRADIANCE PHASE 
TRANSITION 

We now demonstrate that the system considered here 
is capable of producing phonon superradiance. The BEC 
of N two- level atoms with level spacing Hujl can be de- 
scribed by a collective spin S — N/2. In the weak 
atom-phonon coupling limit, the counter- rotating terms 
a^S + and aS- in the usual Dicke Hamiltonian arc usu- 
ally neglected (rotating wave approximation ) (4l| . In 
the strong atom-phonon coupling regime and including 
atom-atom interactions, one can show that the Hamilto- 
nian of Eqn.(8) can be rewritten as, 



uqS x = -2^ X S Z S X + 4Aai5" z , (26) 



{2 1i S z +lo )S v ^0, (27) 

where ujo — W£ +70+71. From Eqn. (|2"Tl) . we observe that 
there are two classes of solutions depending on whether 
S y = or S x = —^0/271. S y = is the usual superradi- 
ant phase in the Dicke model. For S y = 0, we obtain the 
steady state population difference, 



H = Hui r a'a + %{uil +7o)<$ z + hrf%S z 
+ h\(a) +a)(S+ + S-), 



(21) 



where S z , S- and S+ are the collective spin operators 
of the BEC [H. S± = S x ± iS v and S x = V2g F F x 
[38 

A : 

7i 

fvj\Sl {Si ^ 1) is expected to show rich collective dy- 
namics and phase diagrams due to the fact that 71 can 
be both positive and negative. We now study the dy- 
namics arising out of the Hamiltonian of Eqn. ([2"Tj) . The 
equations of motion for a, S z and S- are given by 



(F x is the x component of the atomic spin). Also, 
: G N /V2S, G N = GVN, 70 = (K ee - K gg )/A and 
[K gg + K ee — 2K eg )/2. Note the nonlinear term 



-ioj r a — Kfl 



(22) 



S z 



8A 2 w r -27 1 (K 2 +w2)- 



(28) 



The critical coupling strength for the onset of superra- 
diance is obtained by setting S = (0,0,±7V/2). One ob- 
tains, 



A 2 iV ± =± K±7i^V)(« 2 +^) 



4uj r 



(29) 



/ 2 1 2\ 
_ j_ U)L (U r + K ) 



For 7! = 0, 70 = 0, we have A 2 N± 
which is the usual Dicke model critical point (4lJ. For 
71 7^ 0, we find that the two body atom-atom interactions 
turn out to be a convenient and new handle to tune the 
critical point. The second solution S z = — wo/271 leads 
to a — 0, which turns out to be the normal phase. 



S z = A(a f + a)(5- - S+), (23) 

S- = -i(u)L + 7o + Ji)S- - 2ijiS-S z 

+ i2X{a^+a)S z . (24) 

Here we have taken the phonon damping rate T/2 = k. 
A. Mean field solutions 



B. Fluctuations 

We now discuss the instability of the normal phase and 
the inverted phase by considering fluctuations of phonon 
number and the spin S- around the steady state. To this 
end, we write a — > a + 5a and S- = S— + where 
a = 0, S 1 - = and S z = =pN/2. We consequently obtain 
the linearized equations 

Sa = — (k + iu) r ) Sa — i\(5S- + 5S+), (30) 



We first study the long time, mean-field steady state 
solutions with 5 = and a — 0. Following [42j], we iden- 
tify two stable states; the normal state (J|, denoted by 
N in the phase diagrams), with all spins pointing down, 
S z = —N/2 and no phonons, a = 0, and the inverted 
state (f|-, denoted by I in the phase diagrams) with all 
spins pointing up, S z — N/2 and no phonons. 

We now look for other interesting configurations by 
analyzing the steady states. Writing a = a\ + iai ,S± — 
S x ± iS y , a — 0, S z — and S_ — leads to 

[k + iu! r ]a = —i\S x , (25) 



SS- = -iu)o T SS- T i\N(5a + 60?), (31) 

where ujq t = ljq =p JiN. We look for solutions of the 
form 8a = Ae^ + B* e"!** and SS- = Ce""'* + DV'** 
and equating coefficients with same time dependence, one 
obtains equations for A, B, C and D. Following (42J, 77 
satisfies, 

(ry 2 - u? r - K 2 )(ry 2 - cJ |) T 4A 2 7Va; r a7o T 

- i2nr]{f] 2 -cJ 0t ) = 0. (32) 
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The boundary between unstable (exponentially grow- 
ing ) and stable (exponentially decaying) solutions cor- 
responding to Eqn. ()32j) having real solutions for 77. 
The imaginary part of Eqn. (|32[ ) vanishes when 77 = or 
77 = cJo T . 77 = ujq t implies that the real part of Eqn. (|32"|) 
vanishes when ujq t = 0. This implies both the normal 
and inverted states become unstable when loq — ±71 N. 
For rj — 0, the real part of Eqn. (|32l) becomes zero when 



regions of N phase and / phases increases. 



(uif + k 2 )oJo^ T 4\ 2 Nu} r uio T — 0. This gives the same 
condition as Eqn. (|2T)1) . Consequently this implies that 
the onset of the superradiant phase is accompanied by 
the instability of the normal(iV) and inverted phase (I). 

The dynamical phase diagram corresponding to the dy- 
namics of Eqns. (30,31), can be calculated by the cor- 
responding eigenvalues. In this non-equilibrium setting, 
the eigenvalues are given by, 



±2X 2 N 
Wt = — ± 



4A 4 7V 2 



w T7i^ V( Wo= F7i-^) 2 



(33) 



C. Phase diagrams 

The dynamical phase diagrams emerging from the 
eigenvalues of Eqn. (33) are shown in Fig. 3. For 71 iV = 
and for the positive eigenvalue w_ , the phase diagram of 
Fig. 3(a) reflects the equilibrium phase diagram of the 
Dicke model, having a transition from the normal phase 
(N) at low Xy/N to the superradiant normal phase (SRN) 
at higher value of X^/N. This dynamical phase transi- 



tion occurs at Xy/N 



As LJ- 



0, the 



critical value of Av N required for superradiance tends 
to infinity. Also we find as in (42[, for negative eigen- 
values CcL|_, this open dynamical system shows signature 
of non-equilibrium dynamics; the normal state (Jj.) be- 
comes unstable and an inverted state (S z — N/2, a = 0, 
denoted by /) and superradiant inverted phase (denoted 
by SRI) emerges which is a stable state. The inverted 
state is the mirror image of the normal state which is also 
reflected in the equation of motion Eqns.(22)-(24), which 
have an inversion symmetry for uj r — > — uj r , S — > —S, 
a — » a* and 71 = 0. 

In the presence of a finite 71 , this inversion symmetry is 
immediately broken which is evident from the Eqns. (22)- 
(24). Fig. 3(b) and Fig. 3(c) illustrates this broken in- 
version symmetry. For 71 AT = —0.04 (Fig. 3(b)), the 
phase boundary between the N phase and SRN phase 
recedes to higher Xy/N values while the phase bound- 
ary connecting the I phase and the SRI phase shifts 
towards lower Xy/N values. Exactly the opposite is true 
for 71 N — 0.04 (Fig. 3(c)) A further increase in the val- 
ues of 71 TV = ±0.06 gives rise to regions where both the 
SRN and SRI phases coexist (Fig.3(d) and Fig.3(e)). 
The influence of increasing the phonon damping rate (k) 
is shown in Fig. 3(f). An increase in the phonon damping 
rate separates the SRN and SRI phases further i.e the 
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Figure 3: Dynamical phase diagrams depicting the various 
phases (normal phase N, inverted phase I, superradiant nor- 
mal phase SRN and superradiant inverted phase SRI). Plot 
(a): 7i7V = 0,k = 8.1 , Plot (b): j x N = -0.04, Plot (c): 
7iiV = 0.04, Plot (d): 7i7V = 0.06, Plot (e): 7l N = -0.06 , 
Plot (f): 71 AT = 0, k = 16.1 



IV. CONCLUSIONS 

In conclusion, we have shown that a phonon laser can 
be fabricated based on magnetically coupling a ultra-cold 
atomic cloud to the mechanical oscillations of a nanoscale 
magnetic cantilever in close analogy to a two-level optical 
laser system. By controlling the number of atoms, one 
can switch between solitary pulses and transient pulses. 
The transients can also be controlled by the atomic two 
body interaction. We have also demonstrated that the 
system considered here is capable of producing phonon 
superradiance. For large atom-mechanical mode cou- 
pling, the system can be described by the Dicke type 
Hamiltonian. The two-body interaction gives rise to a 
nonlinear term proportional to S 2 , which gives rise to 
some interesting phase diagrams. By appropriately tun- 
ing the two-body interaction, we get regions in the phase 
diagram where the superradiant normal and superradiant 
inverted phase coexist. 
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